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1. Introduction 


Quantum Chromodynamics (QCD) admits two interesting limits, which can be treated 
by similar methods. In the sector of the heavy quarks (c and b) one observes to 
leading order the so-called heavy quark and the related spin symmetry [1]. The heavy 
quarks can be treated in a non-relativistic framework which is called Heavy Quark 
Effective Field Theory (HQEFT), reminiscent of the well-known Foldy-Wouthuysen 
transformation for heavy Dirac particles. The QCD Lagrangian for the heavy flavors, 
collectively denoted by the held Q, takes the simple form 


— Q{i P —Mq)Q , ( 1 . 1 ) 

where all the propagation of the heavy quarks and their interactions can be handled as 
power expansions in 1/Mq. Similarly, for the light quark sector, one observes that the 
current quark masses are small compared to the typical scale of the strong interactions. 
This dehnes the chiral limit of QCD, which can be analysed making use of chiral 
perturbation theory (CHPT) [2]. This is the effective held theory of the Standard 
Model at low energies, with the pseudo-Goldstone bosons n, K, rj the pertinent degrees 
of freedom. In the presence of matter helds, like e.g. baryons, a complication arises 
due to the large baryon mass scale, which is comparable to the scale of chiral symmetry 
breaking, ms ~ ~ 1 GeV. Resorting to methods taken from HQEFT, Jenkins and 

Manohar [3] showed how to move the troublesome baryon mass term in a string oil /tub 
suppressed interaction vertices. This allows for a consistent power counting since the 
baryon propagator to leading order is given by 

S{u) = ^ {e - 0+) , (1.2) 

with uj = v-kin terms of the four-velocity vector and the small residual momentum 

k, V ■ k rriB- As wanted, the baryon mass has disappeared from the propagator. 

The simplest way to systematically calculate the I/tub corrections was spelled out in 

ref. [4], in which the path integral formalism for HQEFT developed by Mannel et al. [5] 
was extended to baryon CHPT. 

Clearly, in heavy baryon CHPT as well as in HQEFT, the massive degrees of free¬ 
dom behave essentially non-relativistically and it is thus not obvious how to extend 
the notion of wave function renormalization to such a situation. In relativistic baryon 
CHPT, this is not an issue since one can apply standard quantum field theoretical meth¬ 
ods, as detailed in [2] [6]. In the work of ref. [4], the wave function renormalization was 
dehned via the derivative of the nucleon self-energy at cu = 0, leading to a momentum 
independent result for the heavy nucleon Z-factor. Since the propagator eq.(1.2) 
develops a pole at this value of cu, it is the natural point to expand around, but one 
could equally chose other values of u to dehne the Z-factor.A somewhat different 


are much indebted to Thomas Hemmert for clarification on this topic. 



interpretation was given in ref. [7]. A more detailed analysis of this particular aspect 
was performed by Ecker and Mojzis [8], who argued that the Z-factor can not be a 
constant but rather depends (in momentum space) on the chosen frame via the baryon 
momentum. They for the hrst time stressed the role of the heavy fermionic sources and 
within their scheme, the contribution from these sources is entirely given by and 
one thus does not have to perform any explicit calculation for terms involving these 
heavy sources (once the Z-factor is determined). Note that the Z-factor given in that 
paper for the “BKKM” approach is not correct, it should be momentum-independent. 
Such an observation was independently made in ref.[9].^^ This momentum dependence 
is, however, also present in the treatment a la ref. [4]. In that approach, the tree graphs 
are calculated from the relativistic tree level Lagrangian and then expanded in inverse 
powers of the nucleon mass up to the needed accuracy. More precisely, one has to 
include all relativistic tree level Lagrangian terms + ... which, when ex¬ 

panded, can contribute to the order of M.„/mN one is after. This procedure contains 
automatically the momentum dependent pieces through the spinor normalization, as 
explicitely shown for the case of elastic pion-nucleon scattering in [10]. A general proof 
that this method always leads to the correct results has not been given. The aim of 
this paper is to set up a very simple scheme for wave function renormalization in heavy 
fermion effective held theories which parallels as closely as possible the conventional 
quantum held theory approach. As we will show, it is very useful to elucidate the in¬ 
terrelationship between the various approaches found in the literature. It also provides 
us with the proof that expanding the relativistic tree graphs indeed leads to the correct 
result, as hrst conjectured in [11]. The method developed should also be of interest for 
HQEFT, since to our knowledge this issue has not been addressed in detail there. An 
exception to this is ref. [12] in which the question of how to properly normalize spinors 
in HQEFT is discussed. 

The manuscript is organized as follows. In sec. 2 we review the salient features of 
heavy baryon chiral perturbation theory (HBCHPT) necessary to keep our presentation 
self-contained. In sec. 3 we establish a novel scheme to dehne wave function renormali¬ 
zation in heavy fermion EFTs, based on a set of four simple conditions and interpreting 
the so-called light components of the heavy fermions as Dirac spinors. This scheme is 
free of a set of ambiguities, which naturally arise in heavy fermion effective held theo¬ 
ries. Section 4 is devoted to an alternative interpretation, in that the light components 
are treated as two-component Pauli spinors. This allows us to establish the relation 
between the general dehniton given before and the method of expanding the relativistic 
tree graphs, as commonly done [13]. Since S-matrix elements and transition currents 
do, of course, not depend on the way one dehnes wave function renormalization, we 
consider in detail the electric Sachs form factor in sec. 5. We show that all different cal- 
culational schemes lead to a constant contribution from the Born terms in conjunction 
with the appropriate wave function renormalization. We also show (in an appendix) 


are grateful to Gerhard Ecker for confirmation on this statement. 



that the loop graphs do not renormalize the charge, as it should be. The final com¬ 
ments and summary are given in sec. 6. Some technicalities and additional remarks are 
relegated to the appendices. 


2. Heavy nucleon effective field theory 


We briefiy review the path-integral formulation of the chiral effective pion-nucleon 
system. This follows largely the original work of [4], which was reviewed in [13]. The 
interactions of the pions with the nucleons are severely constrained by chiral symmetry. 
The generating functional for Green functions of quark currents between single nucleon 
states, Z[j, ? 7 ,f/], is defined via 


exp =M 


exp i 


Sn-K + riN + Nn ) 


with and denoting the pion and the pion-nucleon effective action, respectively, 
to be discussed below, rj and fj are fermionic sources coupled to the baryons and j 
collectively denotes the external fields of vector (n^), axial-vector (a^), scalar (s) and 
pseudoscalar (p) type.^^ These are coupled in the standard chiral invariant manner. 
In particular, the scalar source contains the quark mass matrix M., s{x) = M. + .... 
The underlying effective Lagrangian can be decomposed into a purely mesonic (tttt) 
and a pion-nucleon (ttIV) part as follows (we only consider processes with exactly one 
nucleon in the initial and one in the final state) 


= d^-K-K + 


ttN 


( 2 , 2 ) 


subject to the following low-energy expansions 


£^,, = £( 2 )+£( 4 ) + .. £^^ = £«+£ 


■'irN 


,( 2 ) 

■'nN 


C 


( 3 ) 

irN 


+ 


wN 


(2.3) 


where the superscript denotes the chiral dimension. The pseudoscalar Goldstone fields, 
i.e. the pions, are collected in the 2x2 unimodular, unitary matrix U{x), U{(j)) = 
= exp{i(j)/F} with F the pion decay constant (in the chiral limit). The external 
fields appear in the following chiral invariant combinations: , 

and X = 2Bo (s ip). Here, Bq is related to the quark condensate in the chiral limit, 
Bq = K0|?5'|0)|/F^. We adhere to the standard chiral counting, i.e. s and p are 
counted as 0{q‘^), with q denoting a small momentum or meson mass. The effective 
meson-baryon Lagrangian starts with terms of dimension one, 

O 

= + , (2.4) 

^^The external vector field should not be confused with the four-velocity to be dehned later on, 
which is denoted by the same symbol. 


with mo the nucleon mass in the chiral limit and = i[u^{d^ — ir^)u — u{d^ — 

The nucleons, i.e. the proton and the neutron, are collected in the iso-doublet T, 
= (p,n). Under SU{2)l x SU{2)r, T transforms as any matter held. denotes 
the covariant derivative, + T^T and T^ is the chiral connection, T^ = 

I [u\d^ — ir^)u + u{d^ — il^)u^]. Note that the hrst term in Eq.(2.4) is of dimension one 
since {ip—mo) T = 0{q) [14], The lowest order pion-nucleon Lagrangian contains two 

O 

parameters, namely mo and 9a- Treating the nucleons as relativistic spin-1/2 helds, 
the chiral power counting is considerably complicated due to the large mass scale mo, 
(9o 4/ ~ mo 4/ ~ A^4/, with ~ 1 GeV the scale of chiral symmetry breaking. A 
detailed analysis of this topic can be found in [6]. This problem can be overcome in 
the heavy mass formalism proposed in [3]. We follow here the path integral approach 
developed in [4]. Dehning velocity-dependent spin-1/2 helds by a particular choice of 
Lorentz frame and decomposing the helds into their velocity eigenstates (also called 
’light’ and ’heavy’ helds), 

Hy{x) = exp{imon ■ x} Pp ^(a;), hy{x) = exp{imon ■ x} Pp 4/(a;) , (2.5) 

the mass dependence is shuffled from the fermion propagator into a string of 1/mo 
suppressed interaction vertices. The projection operators appearing in Eq.(2.5) are 
given by 

P± = i^, P:H = H,P;h = h, P+ + F- = l, (2.6) 

with the four-velocity subject to the constraint = 1. To be specihc, the nucleon 
four-momentum has the form 

= moV^, + kf, , (2.7) 

where is a small residual momentum, v ■ k mo. In the basis of the velocity 
projected light and heavy helds, the ehective pion-nucleon action takes the form 

SttN = J d'^x AP[jj-hyC hy+ hy B p[y + Hy'joB^'Johj^ . (2.8) 

The matrices A, B and C admit low energy expansions, e.g. 

A = A^^'^ + + A^^) + ... , (2.9) 

and similarly for B and C. Explicit expressions for the various contributions can 
be found in [13]. Furthermore, we split the baryon source helds ri{x) into velocity 
eigenstates, 

Ry{x) = exp{imov ■ x} Pp ri{x ), Pv{x) = exp{imov ■ x} Pp ri{x) , (2.10) 

and shift variables, hy = hy — C~^ {B Ply + p„), so that the generating functional takes 
the form 

exp[iZ] = A/”A/i j[dU][dHy][dHy]ex-p{iST,T, + iS'.^j^} (2.11) 



in terms of the new pion-nucleon action S'^^, 

SIj, = Jd^x!^ HM + 7o5t7o C-^B )H, + H,{R, + ^,B^^oC-^Pv) 

-\-{Rv PvC ^ B)H^-\-p^C ^Pt,| . (2.12) 

The determinant A/^ related to the ’heavy’ components is identical to one, i.e. the pos¬ 
itive and negative velocity sectors are completely separated. The generating functional 
is thus entirely expressed in terms of the Goldstone bosons and the ’light’ components 
of the spin-1/2 fields. The action is, however, highly non-local due to the appearance 
of the inverse of the matrix C. To render it local, one now expands C~^ in powers of 
1/mo, he. in terms of increasing chiral dimension. To any hnite power in 1/mo, one 
can now perform the integration of the ’light’ baryon held components by again 
completing the square, 

K = H, - T-^ {R, + ^oB^ 'yoC-^ Pv), T = A + 70 70 G-'5 . (2.13) 

Notice that the second term in the expression for T only starts to contribute at chiral 
dimension two. Finally, we arrive at 

exp[iZ] = A/”' J [dU] exp{iS'^ 7 r + } , (2-14) 

with A/”' an irrelevant normalization constant. The generating functional has thus been 
reduced to the purely mesonic functional. is given by 

ZnN = -J d^x I p„ (G-i B T-i 70 7oC'-i _ ^-1 ^ 

+ p„ ( C-^ B T-i )R, + R, (T-i 70 B^ 70 ) p. 

+ R,T-^R, I . (2.15) 

At this point, some remarks are in order. First, physical matrix elements are always 
obtained by differentiating the generating functional with respect to the sources rj and 
fj. The separation into the velocity eigenstates is given by the projection operators as 
dehned above. As shown in ref. [ 8 ], the chiral dimension of the ‘heavy’ source p^ Pv V 
is larger by one order than the chiral dimension of the ‘light’ source, R^ ~ P^V- The 
effective Lagrangian can be readily deduced from this action. For later use, we give the 
hrst two terms, following the dehnitions of [13], 

4v = H,[iv -D+h S-u^H, , 

o2 

+ (C2 - (V ■ !l)^ + C3'U ■ M + . . , 


( 2 . 16 ) 



where the ellipsis stand for three other terms not needed for the following discussions. 

is the covariant spin-operator a la Pauli-Lubanski, 5*^ = | 75 v^, subject to the 

constraint 5* ■ n = 0 and we work in the isospin limit rUu = rrid- Traces in flavor space 
are denoted by Notice that the spin-matrices appearing in the operators have all 
to be taken in the appropriate order. The explicit symmetry breaking is encoded in 
the matrices x± = ± The q are hnite low-energy constants (LECs), their 

values have been determined [15] [16] [10]. 

At one loop level, divergences appear. These can be extracted either by direct 
Feynman graph calculations or, more elegantly, directly from the irreducible generating 
functional [17] [18]. Here we give the form relevant to fourth order in SU(2), the details 
can be found in [19] 


Zirrli,-RJ = Jd'‘xd‘x'd‘yd'‘fi'R„{x) !/) - y') 

+ £<'>+'">(!,, y') + T,f\y,y') ] x’) R.(x') (2.17) 

in terms of the self-energy functionals Si 2 , 3 , and denotes the classical propagator 
in the presence of external helds. Si refers to the self-energy graphs at order and 
the same diagram with one dimension two insertion on the nucleon line. S 2 collects 
the tadpoles at orders and and S 3 refers to the dimension two vertex corrected 
self-energy diagrams. The whole machinery and the complete fourth order counterterm 
Lagrangian is spelled out in [19]. 

So far, there exist three different approaches how to calculate matrix elements. 
The hrst one is based on ref. [4] and will be referred to as “BKKM” in what follows. It 
amounts to a “hybrid” calculation. The tree graphs are worked out from the relativistic 
pion-nucleon Lagrangian and then expanded in powers of l/rriN to the order one is 
interested in. For the reasons mentioned above, the loop graphs are calculated in the 
heavy nucleon framework. In particular, the light helds H are treated as Pauli spinors 
and the corresponding Z-factor is entirely given by the loop graphs and is momentum- 
independent. This method is very convenient for calculations and gives the correct 
results to orders q^ and as will be shown later. The disadvantages of the method 
are twofold. First, such a hybrid approach does not appeal to everybody and second, 
it is not clear how it can be extended correctly to higher orders. Second, Ecker and 
Mojzis [ 8 ] have set up a scheme which stays entirely within the heavy fermion approach, 
however, matrix elements are matched to the corresponding relativistic ones. This 
method should be applicable at any order. The derivation rests on the interpretation 
of the H helds as Dirac spinors. The tree graphs are all calculated within HBCHPT. 
It is important to note that a different Lagrangian is used as compared to BKKM. 
This Lagrangian is subject to held transformations to eliminate the equation of motion 
terms and can be found in [20]. In this framework and for the Lagrangian of [20], the 
wave function renormalization is momentum-dependent (with the exception of forward 



matrix elements), 


4a3M2 

Zn{Q) — 1 h- o —I" 




m% 


Ami 


+ ... 


(2.18) 


N 


where the residual momentum Q is dehned in terms of the on-shell nucleon momentum 
Pat and its physical mass = mp^v + Q. The LEG 03 is related to Ci in eq.(2.16) 

via 03 = mpj Cl- This method has the advantage of staying within one given held 
theoretical framework, however, the Z-factor they give for the BKKM approach should 
be momentum-independent. For a more detailed discussion of this approach and the 
matching to relativistic matrix elements, we refer to ref. [ 8 ]. Third, a variant of the 
BKKM approach, which can easily be extended to higher orders, has been proposed 
by Fettes et ah [10] (called FMS from here on). Again, the light helds are treated as 
Pauli-spinors. The tree (Born) graphs are, however, calculated in the heavy baryon 
limit and the Z-factor consists of two pieces. 






(2.19) 


with M the relativistic spinor normalization. 


M = 


rriN 


2m 


N 


( 2 . 20 ) 


where Ep = yp'^ + m'j^ denotes the full relativistic nucleon energy. For the case of pion- 
nucleon scattering to order 0{q^), it was demonstrated that this method reproduces the 
result from expanding the relativistic tree graphs, i.e. the BKKM approach. In fact, 
as we will show later, one can show this equivalence quite generally for all one-loop 
processes including fourth order. This method has the advantage that it also stays 
within HBCHPT and is thus not a hybrid type of calculation. However, the important 
normalization factor A/"^ which enters the Z-factor, eq.(2.19), is not directly given by 
the heavy baryon theory. It is important to note that in the rest-frame Vp = (1,0) a 
connection between the relativistic and the Pauli-spinor interpretation is given through 
the relation 


P,tu 



( 2 . 21 ) 


with u (x) denoting the conventional four(two)-component Dirac (Pauli) spinor. All 
three methods used so far share one further common disadvantage. One can only 
calculate the Z-factor if one approaches the physical pole in the direction of the four- 
velocity V, i.e. when the nucleon four-momentum p defined in terms of the nucleon 
mass in the chiral limit, mo, is taken to its on-shell value p^ dehned in terms of the 
physical mass mpp, the difference is {pn — p)p ~ Vp. Clearly, it would be preferable 
to have a dehnition of the wave function renormalization that does not depend on 
how one approaches the physical on-shell point. In particular, in the Pauli-spinor 



interpretation one selects the rest-frame as the preferred frame and that makes the 
Z-factor depending on this particular choice of the four-velocity vector. The aim of 
the present work is to set up a scheme which does not have the various short-comings 
discussed so far and allows to give a very simple and concise dehnition of wave function 
renormalization following as closely as possible conventional quantum held theories. 
Furthermore, as we will demonstrate, in this novel scheme it is particularly easy to 
compare the existing approaches and pin down the pertinent differences. Clearly, since 
the Z-factor is not an observable, one is free to chose ones own dehnition, in particular 
also the momentum, about which one expands. The only condition to be fulhlled is to 
work consistently within the chosen framework. 

3. Wave function renormalization reconsidered 

In this section, we wish to establish a scheme, allowing for a dehnition of the heavy 
fermion Z-factor, which we call from here on, subject to the following four condi¬ 
tions: 

1) The dehnition of should be independent of the choice of the four-velocity 
vector V. 

2) Its dehnition should only involve the physical helds. 

3) At tree level, one should have = 1. 

4) The dehnition of should be independent of the way one approaches the phys¬ 
ical on-shell momentum, p ^ Pn- 

As for the last point, we remark that this is also stressed in ref. [8], however, for the 
actual calculation a particular “clever” choice was made to control the on-shell limit. 
While that is certainly legitimate, our aim is to avoid such a choice from the beginning. 
Although we will establish this scheme for the particular case of two havor heavy 
baryon chiral perturbation theory, the method is more general and can be applied to 
other situations. Of course, we stress again that Zj^ is not an observable and that 
many alternative schemes exist to dehne it. As we will argue, the one presented here is 
particularly useful to shed light on the various calculational approaches used so far and 
helps to clarify the pertinent interrelationships. A precise dehnition of Ztv is relegated 
to app. A. 

The starting point of our discussion is the interpretation of the light helds H as Dirac 
spinors, following ref. [8]. In this case, we have to consider all helds in the generating 
functional. Thus, instead of splitting the sources in light and heavy components, it is 
advantageous to work with the physical external sources. We thus slightly rewrite the 
generating functional, eq.(2.15), as (the projection operators are kept for clarity) 

Z = - j d^xfi { P+(A + B'C-^B)-^P+ + P-C-^B{A + B'C-^B)-^P+ 


o 



+ P+(A + B'C-^B)-^B'C-^P- (3.1) 

+ P-C-^B{A + B'C-^B)-^B'C-^P; - P-C-^P-] r) , 

with 

(3.2) 

and the ’tilde’ reminds us that the factor expjimon ■ x} is still included in the external 
sources, r) = exp{imon ■ x} rj. The physics is given by the Green functions, i.e. the 
derivatives with respect to the external sources. For S-matrix elements, we only need 
to concern ourselves with the poles of the Green functions when the particles are on 
mass shell. In the generating functional, eq.(3.1), only the term [A + B'C~^B)~^ can 
give rise to a pole as long as one considers small external momenta. To specify the 
most general n-point Green function, we dehne 

T := A + B’C-^B = To + Tj , (3.3) 

where the index ’0’ indicates the absence of external helds and the subsript ’/’ denotes 
the interaction matrix, which to lowest order is nothing but the pion coupling QaS-u. 
The inverse of T is given by 

r-‘ = (r„ + T,)-' = ry(i + 

= r„-‘ - - G„ . (3.4) 

Note that T is nothing but the amputated amplitude for a general process. Stated 
differently, it amounts to all Born graphs with m interaction vertices and m — 1 inter¬ 
mediate nucleon propagators. The pertinent n-point Green function Gn thus takes the 
form 

G„ = P+ G„ P* + P- G„-‘B„ G„ P+ + P+ G„ P;G„-‘ P" 

+p-c„-‘b„g„b;c„-‘p- 

= (i + p+ G„ p+ (i + s;c„-‘) 

= (i + G„-‘B„) r„-‘ P+ fp+ (i + b;g„-‘) . ( 3 . 5 ) 

To arrive at this result, we have used the commutation relations between the operators 
A, B and C and the projection operators P^. We note that eq.(3.5) agrees with eq.(ll) 
of [8]. The S-matrix is now given by reinstating the external legs, 

5 = h(pO(/ - m^)(l + Go-'Po)To-'P+TP+To-'( 1 + P'Go-')(^ - m^)u{p) . (3.6) 

Let us hrst calculate Tq”^, 

^0 = ^0 + BqCq ^ Bq 

= vk + ^^(2mo + vk- + dciM^ 

= (2mo + vk — 4ciM^)“^ (2movk + k'^ + SmoCiM^ — 

= (2mo + n/c — 4ciM^)“^ — ttIq + SmoCiM^ — , (3.7) 



where /c^ is the small residual momentum and — [v ■ k)v^ is the momentum 

orthogonal to the direction given by the four-velocity v. is the leading term in 
the quark mass expansion of the pion mass, = M^[l + (9(g^)]. One can factor out 
the term Cq^ since it does not have any pole as long as one restricts oneself to small 
momenta. For a given hxed order, the corrections obtained by this factorization are 
always of higher order and can thus be neglected. Conseqnently, only the second term 
in eq.(3.7) contains a pole, which leads to the well-known mass shift 

6m = —4ciM^ . (3.8) 

So we are left with the calculation of 

- mN)il + Cq^Bq)Tq^ = - mN)iCo + Bq){p‘^ - 

= - mN){^ + mN + vk{l- ]6)){p^-m%y 

= l + 2{^-mN)vky-m%yP~ . (3.9) 

Using furthermore = 0, the S-matrix follows as (see also app. A) 

5 = u{p')pypy{p) = h(p')^vP+TP+zjvu(p) 

= h(p')phys \[^Pv T P+ y^M(p)phys , (3.10) 

which means that in this case the Z-factor is exactly one, = 1. Of course, there are 
still corrections from the loops, which will be evalnated subseqnently. We stress here 
that this interpretation allows for a clear and concise dehnition of the Z-factor in that 
only the loop graphs lead to a non-trivial contribntion. 

We now consider the effects of pion loops. For that, we expand aronnd the classical 
solution of the fermion propagator in terms of pionic fluctnations [6]. This means for 
the matrices A,B and C dehned in eq.(2.8), 

A ^ + A^^) + + ... , 

B ^ B^y ^(1) + ^(2) +..., 

c^y + yyyy..., ( 3 . 11 ) 

where the superscript ’(i)’ counts the nnmber of pionic flnctuations^^ (for details, see 
e.g. refs. [6] [18]). Note that we do not need to consider baryonic fluctuations since we 
only work at small momenta. To calculate the corresponding mass shift, we only need 
to work ont the inflnence of these flnctnations at the pole, 

r„-‘ = K> + r„'‘> + rf +.. .]-■ 

= +... , ( 3 . 12 ) 

^"‘These should not be confused with the chiral dimension used before. 



where (5'^^)“^ = Tg*. Integrating over the fluctuation flelds, the second term in eq.(3.12) 
vanishes, whereas the third and the fourth give the tadpole and the self-energy con¬ 
tribution, respectively. We can thus bring the expression for the self-energy into a 
compact form 

Tg-^ = + ... = [Tg^* - S]-1 . (3.13) 

The last equation is exact in the sense that if one were to calculate the contributions 
from all irreducible graphs to S, the reducible ones follow from the geometric series (as 
it is well-known). The inverse propagator follows as 

S~^ = vk + k^{2mo + vk — 4ciM^ — Sb^M^Y^k^ 

+4ciM2 + 86oM^ - t^^\vk) - t^^\vk) 

= (2mo + vk — 4ciM^ — 

X — mg -I- 2mo(4ciM^ -|- SboM"^) — (4ciM^ -|- SboM^Y 

-{2mo + vk-AciM‘^-SboM%Y^\vk) + t^^\vk))] , (3.14) 

with (u) denoting the third and fourth order loop contribution to the nucleons’ 
self-energy. Here, 6o is a combination of the LECs of three counter terms, which 
in fact lead to a quark mass correction of the LEG Ci. The numerical value of &o 
is at present not known. It could be obtained from a fourth order analysis of the 
baryon masses and a-terms (for a model-dependent determination within SU(3) baryon 
CHPT, see e.g. [22]). In what follows, we absorb the contribution of the bo term via 

c[ Ml = ciM^ + 2boM^ . (3.15) 


Explicit calculation using the usual procedure of renormalization in CHPT [2], leads 
to the third order contribution to S, 


= eSp(o.) + Sgi(o.) 
_ ^9A . .2\ f ^ 




arccos ■ 


.29AYMl-2Y)^ln- 

2F2 ^ ^ ’ 16n^ 


+ - SMlud^sY) . (3.16) 


Similarly, the fourth order self-energy reads 
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(3.17) 


with A the scale of dimensional regularization. All quantities have been set to their 
physical values since the differences to the chiral limit values only appear at next order. 
Of course, eq.(3.16) agrees with the third order self-energy expression given in [4]. The 
logarithms in eq.(3.17) stem from the self-energy graphs (fourth last) and the tadpoles 
(third last line), respectively. The 621 199 ^-re renormalized fourth order LECs taken 

from table 1 of ref. [19] (note that the fourth order LECs are called d* in [19]. To avoid 
confusion with the labelling in the FMS Lagrangian, we call them bi here). By a proper 
redehnition of the renormalized fourth order LECs, one can absorb all the logarithmic 
terms, 

i>KA) =-A h + In ■ (3.18) 


Consequently, all ln(Mn-/A) terms vanish in eq.(3.17). The corresponding di and bi are 
all zero, as detailed in [10]. With a similar procedure for the d*, one can also absorb 
all the logarithms in eq.(3.16), for details see [10]. From the pole position we read off 
the mass shift 


5m = — 


2mo + 5m 
2mo + vk — 4c( 
2mo + 5m 


4c' Ml +-4c( Ml 

^ ^ 2mo +5m ^ " 


{t^^\vk) + t^^\vk)) 


(3.19) 


To proceed, we have to work out the value u/c|phys, 

- ml - k‘^ k'^ 


+ ... . 


(3.20) 


Let us comment on the /c^-dependent terms. The third order self-energy does not 
depend on vk, since vk = O^q"^), i.e. such terms can not contribute at third order. At 
fourth order, however, one gets a term of the type uA; S*^^^'(0). The nucleon mass shift 
is, of course, not /c^-dependent. At fourth order, all terms ~ cancel, as we will make 
explicit below. The fourth order mass shift reads 
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with the pertinent /9-functions, {At^FY l3{h\Q^) = —95f^/(16mo) and {At^FY I3{d2d) ~ 
—^g\/lQ. As promised, the k"^ terms add up to zero. In that basis, the mass shift to 
fourth order takes the form 


5m = -Ac[Ml 


^g\Ml 3Mtc, 
327rF2 ^ 1287r2F2 




+ 0{qY . 


(3.22) 


Here, all masses and couplings are set to their physical values, the error made by this 
procedure is of higher order. It is instructive to work out these corrections numerically. 
Using the LECs as determined in [15] (and ignoring the quark mass renormalization of 
Cl, i.e. setting c( = ci), we get 


= (72.5 - 15.1 + 0.4 - 0.4) MeV , (3.23) 


which shows that with the exception of the (undetermined) (hidden in c(), the 

fourth order corrections are tiny. The general structure of the fourth order contribution 
to the nucleon mass was already given by Kallen [21], but we do not agree with some 
of her coefficients. 

We now come back to the Z-factor. We have to generalize eq.(3.5) in that all 
possible loop effects have to be taken into account. In a symbolic language, this reads 

G„ = (1 p:fp: T„-‘(1 + b;c„-‘) 

+ (1 + P+fP+V(1 + BJC„-') 

+ (1 + pY f PY T-\1 + 

' -V-' 

+ (1 + p*fp: T„-'(1 + 

'-V-" 

+ ... (3.24) 


where the ’underbrace’ indicates that loops have to be constructed from the pertinent 
operator structures involving Cq^ and Bq. A closer look at these contributions reveals 
that only the hrst two underbraced combinations of operators have the correct pole 
structure, however, they only start to contribute at order (when one counts the 
small momenta at the level of the Lagrangian). This means that loops involving the 
heavy sources do indeed only start to contribute at that order. Consequently, for 
calculating the Z-factor to 0{qY we have to work out (in analogy to eq.(3.9)) 

{Y - m/v)(l + C'o-'Ho)[To - S(o^)]-i 
= (^-m/v)(C'o + Ho)[C'oTo-C'oS(o^)]-i 

= {Y ~ ^n){Y F ^0 — + 2vkP~)[CoTQ — CqT^{uj)]~^ , (3.25) 

with PYPF — 0- is the value at the (physical) pole, 
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where the prime denotes differentiation with respect to p and we used 



rriN 

mo 
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Finally, we arrive at 
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Note that the hrst correction from two loops appears at order q^ since the Z-factor is 
the derivative of the self-energy and thus one power of small momentum is absorbed 
in the dehnition of Zn- We remark that the Z-factor is momentum-independent and 
that, by construction, it fulhlls all four requirements spelled out in the beginning of 
this section. We are now in the position to analyze the various methods used so far 
in explicit HBCHPT calculations, mostly based on a Pauli-spinor interpretation of the 
light helds. To end this section, we note that in ref. [8], the Dirac interpretation was also 
used. In contrast to what is done here, held redehnitions have been used to bring the 
effective Lagrangian into a minimal form. We remark that the question of the inhuence 
of such held redehnitions needs to be addressed in more detail. Such an investigation 
is underway but goes beyond the scope of the present paper. 


4. Pauli spinor interpretation 

To clarify the approach used by BKKM and also by FMS, we now interpret the light 
helds H as two-component Pauli spinors. Clearly, this dihers from the interpretation 
given in the preceeding section and the crucial point will be the discussion of the 
contribution from the heavy helds/sources, which in BKKM are hidden in the expansion 
of the relativistic tree graphs and in FMS come in via the spinor normalization factor. 

To be more precise, we reiterate the statement already made in sec. 2 that in 
the rest-frame = (1,0) a connection between the relativistic and the Pauli-spinor 
interpretation is given through the relation u = A/" (y, 0)^, with u (y) denoting 
the conventional four(two)-component Dirac (Pauli) spinor. In light of these remarks. 



consider again the generating functional. In the Pauli spinor interpreation, the pion 
nucleon action has no heavy sources, cf. eq.(2.12), 


S'^j^ = I d^x H,{A + B' C-^B )H, + H,R, + R, H, . 

Completing the square is achieved by setting 

H', = - T-^ R, , T = A + B'C-^ B , 


(4.1) 


(4.2) 


leading to 

Z = - j d^xR^{A + B'C-^B)-^R^ , (4.3) 

which means that the other components do not play any role. The matrices A,B and C 
are the standard ones when the relativistic Lagrangian is turned into its heavy fermion 
form as explained in section 2. Consequently, the inverse propagator is given by 

= Ao + B'^C^^Bo = To , (4.4) 


with To given in eq.(3.7). Therefore, to this order one again obtains the standard mass 
shift 5m = —4ciM^. The loop corrections to this result will be discussed below. It is 
instructive to hrst consider the Z-factor, which is dehned as follows. 
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(4.5) 


where “phys” means that the expression has to be evaluated at the physical value of 
vk. Setting now v = (1,0, 0,0) (i.e. considering the rest-frame of the heavy nucleon) 
gives 

m]\f + Ep tA a\ 

Zn = — -= M , (4.6) 

with M the normalization factor of the relativistic spinors, already given in eq.(2.20). 
This underlines the observation made by Fettes et al. [10], namely that retaining the 
spinor normalization allows to recover the expanded full relativistic tree result directly 
from the heavy nucleon approach (in that paper, pion-nucleon scattering was inves¬ 
tigated). Here, this observation is clearly more general. In fact, the dehnition of the 
Z-factor given here fully justihes the method used by BKKM (at least to one loop order 
and Higher order calculations have not yet been attempted, with the exception 
of [22] [23] and [9]). 

The self-energy calculation including the one loop effects proceeds as outlined be¬ 
fore. We only want to stress that in the Pauli spinor interpretation one has the exact 



same pole and thus the same mass shift as in the Dirac case discussed in the previous 
section. The Z-factor is, however, different. 
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Comparison with eq.(3.29) gives 
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(4.10) 


This shows that the factorization, which appeared somewhat ad hoc in eq.(2.19), is 
exactly reproduced. It has its origin in the precise matching of the heavy baryon to 
the fully relativistic theory. This is mirrored in the momentum-dependence of 
in eq.(4.9) and it justihes a posteriori the methods employed by BKKM (at one loop 
order) and FMS. A precise statement about what happens at higher orders can only 
be made when one performs explicit calculations. This is not the aim of this paper. 

Finally, we perform a calculation based on the Lagrangian given in [20] in the Pauli 
spinor interpretation (we stress that this is not what was done in [8]). Field redehnitions 
have been used to bring the Lagrangian in a minimal form. This naturally changes the 
propagator, to third order it reads 




vk + 


2mo 


mo 
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(4.11) 


with the LEG as related to ci via as = ttin ci. Note that we can not give the fourth 
order result since the corrections due to the held redehnitions have not yet been worked 
out at this order. The mass shift is, of course, identical to the one given so far, but the 
Z-factor is diherent. 




Note that in this formulation Zn is momentum-independent because it has been worked 
out for Pauli spinors. We hnally remark that for HBCHPT, the Pauli interpretation 
can be considered the natural framework. 



5. The charge form factor of the nucleon 


In this section, we explicitely calculate the isovector electric Sachs form factor of the 
nucleon to third order in the chiral expansion. This serves to illustrate the point that 
all methods discussed so far lead to the same result if applied correctly. It also shows in 
a very transparent way the differences of the various frameworks in intermediate steps. 
A detailed discussion of the nucleon form factors can be found in [24]. 

To be specific, consider the nucleon matrix element of the isovector component of 
the quark vector current V]] = ?7^(TY2)g in the Breit frame. It was already shown in 
ref. [4] that in HBCHPT this is the natural frame (since the nucleon essentially behaves 
as brick-wall with respect to the incomimg soft virtual photon). In the rest-frame 

= (1,0), the matrix-element of the isovector-vector current takes the form (note 
that the superscript 'v' refers to the isovector current, not the four-velocity) 


W)iK(o)|]v(p)) = y 


niN 


Xi X 72^71 , (5.1) 


where y is a Pauli spinor with isospin component rj and = {p' — pY < 0^^ is the 
invariant momentum transfer squared. and are the isovector electric 

and magnetic Sachs form factors. We remark that we can replace the Pauli-Lubanski 
spin-vector by the Pauli spin matrices, since = (a, d)/2 is restricted to the two 
upper components. Obviously, we also need this matrix element sandwiched between 
Dirac spinors, in which case it takes the form 


{N{p')\vmN{p))=u{p')p: 


Gl{e)v, + - 

rriN 


p+'u(p) X 1)5— r;i , 

(5.2) 

where u{p) is a Dirac spinor. Here, G'%{k‘^) and G'l^{k‘^) are related to the isovector 
electric and the isovector magnetic Sachs form factor, via 
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We are now in the position to calculate these form factors. Note that we calcu¬ 
late the objects ^ within the Dirac spinor framework whereas the Gg follow in 
the Pauli spinor interpretation. From here on, we concentrate on the charge (elec¬ 
tric) isovector form factor. The pertinent kinematics for the tree level photon-nucleon 
coupling, cf. fig.l, is 


p = mov + Pi, p' = mo V + P2, Pi = E', 


E\ + 


k 


= (^0, , n ■ /c = 0 , V ■ Pi = V ■ P2 = E' = E — mo = 0{q^) 


(5.4) 


where E' denotes the residual energy and pi 2 the soft residual momenta of the in- and 
out-going nucleon, respectively. We can rewrite the isovector Sachs form factors as 

^^The photon four-momentum k should not be confused with the same symbol previously denoting 
the nucleons’ small residual momentum. 



■ (Born terms + loops) , (5.5) 

with the loop contribution being the same in all approaches. 
Their explicit calculation is relegated to the end of this sec¬ 
tion. For the moment, we concentrate on the Born (tree) terms 
and the Z-factor, which are different in all schemes but when 
combined, should lead to the same result, symbolically 

as it follows from the relativistic calculation [6]. Note that from the Born terms we 
only consider the ones with fixed coefficients in the 1/mAr expansion since all others 
contribute to the anomalous magnetic moment, the charge and magnetic radii and so 
on. 

We now calculate explicitely the Born contribution to to third order in small 

momenta within the Pauli spinor interpretation (which is equivalent to FMS approach 
to this order). We first collect all tree level terms free of LECs that can contribute 
at third order to the generic diagram shown in fig. 1. The pertinent operators and 
respective Feynman rules (derived from the FMS Lagrangian) for the photon-nucleon 
coupling are 
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Figure 1: Generic 

tree diagram for the 
photon coupling to the 
nucleon. 
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with the polarization vector of the photon and Q = e{l + t^)/2 = ediag(l, 0) the 
nucleon charge matrix. It is important to note that i-6- ^ calculation 

with a term of dimension g” from the Lagrangian gives at order (since the 
polarization vector is taken out). With these rules, we can straightforwardly calculate 
the Born terms corresponding to the electric Sachs form factor, 
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since pf = P 2 = ~ k /4 and E' is of chiral order two. Expanding now the Born 

Z-factor, see eq.(4.6), to third order in momentum. 


Z 


Pauli 

N 




1 + 


16 m^ 


+ ... 


(5.12) 


we hnally get for the isovector electric Sachs form factor. 


= J\[^ ■ (Bom tcrms) = e (1 + 0{q^)) , (5.13) 


i.e. the tree graphs give a constant form factor with the correct charge to the order one 
is working. We now turn to the BKKM approach. It formally amounts to calculate 

G^(bkkm)('^ 2 ^ = ■ (rel. Born terms + loops) . (5-14) 


As already stated, the calculation of the Born terms is trivial since the relativistic 
Born graphs lead to a constant G^(/c^) and thus for this application, the calculation 
is considerably simpler than using HBCHPT to work out the tree graphs in the Pauli 
spinor interpretation as detailed above. 

Finally, we are left with the EM version. Here, the Born terms take yet another 
form since they start from a different Lagrangian, 
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Multiplying this with the Born term Z-factor as given in eq.(4.12), one hnds that due to 
the term ~ mo/m^ the charge form factor is again constant and properly normalized. 
We add a few remarks: First, the difference between the BKKM and FMS schemes is 
5G^^ ~ (Z^™^ — Z™®) loops ~ 0{q*), i.e. one would probably have to modify the 
BKKM scheme for the wave function renormalization when it comes to two-loop (or 
higher order) calculations. We do not follow this subtlety here in detail, but it should 



be kept in mind. Second, we turn our attention to the Dirac spinor interpretation, 
which allows for comparison of all schemes. The calculation of the Sachs form factors 
can be summarized as follows: 
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which means that the Pauli und the Dirac interpretation give exactly the same result. 

In summary, there are many ways to arrive at the correct result, however, given a 
certain scheme, one strictly has to stay within its rules. It remains to be shown that 
the loops do indeed not renormalize the charge, as it follows from gauge invariance, but 
only lead to a momentum dependence of the electric Sachs form factor. Note that the 
calculation of the loop contribution is the same in all three schemes, provided one has 
a prescription how to treat the spin vector in d dimensions. This calculation was first 
performed in [4] (for the isovector Dirac form factor). In the appendix B, we show a 
different and somewhat unusual way to arrive at the same result. Finally, we remark 
that the considerations presented in this section can easily be extended to fourth order. 
For the sake of brevity, we do not spell out the pertinent details here. 


6. Summary and conclusions 

In this paper, we have considered the questions surrounding wave function renormali¬ 
zation in heavy fermion effective held theories. As an example, we have studied heavy 
baryon chiral perturbation theory for two havors and in the isospin limes. The pertinent 
results of this investigation can be summarized as follows: 

i) The most natural and economic way of dehning wave function renormalization 
in heavy fermion effective held theories rests on the interpretation of the light 
components of the heavy helds (like e.g. the nucleons) as Dirac spinors. In 
that way, one can dehne the Z-factor subject to four conditions as detailed in 
section 3. In particular, the so-dehned Z-factor is momentum-independent and 
its tree graph contribution is equal to one. Furthermore, this prescription can 
be extended to higher orders, i.e. beyond one loop, without problems. Such an 
interpretation is mandated by the correct matching of the heavy fermion EFT to 
its relativistic counterpart. 

ii) All calculations performed so far in HBCHPT have been done under the assump¬ 
tion that the light components of the heavy helds are to be interpreted as Pauli 
spinors. We have shown the equivalence between such an approach and the Dirac 



spinor interpretation, provided one works in the rest-frame v = (1,0) in the Pauli 
case. This allows to justify a posteriori the methods employed by BKKM (l/rriAr 
expansion of the relativistic tree graphs independent of the spinor interpretation) 
and by FMS (inclusion of the explicit four-dimensional spinor normalization in 
the tree graphs calculated from HBCHPT). We also pointed towards some po¬ 
tential complications which arise when one employs held redehnitions. 

hi) When applied correctly, all these different schemes lead to the same physics. As 
an example, we have shown how the tree result for the proton charge form factor, 
= e = const., with e the proton charge and k'^ the photons’ four- 
momentum squared, emerges in the various calculational schemes. We have also 
discussed the non-renormalization of the electric charge due to the loop graphs. 

iv) Furthermore, we have studied the nucleon mass shift to fourth order in the pion 
mass. Apart from an unspecihed counter term contribution, which formally 
amounts to a quark mass correction of a dimension two operator, these corrections 
are tiny. 

We hope that with this paper, the long-standing question of wave function renormali¬ 
zation in heavy fermion effective held theories can dually be put to rest. 
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A. Precise definition of the Z—factor 

We first repeat the standard steps to get from the generating functional Z[j,r],f]] for 
bosonic fields (like e.g. the pions) and fermions (like e.g. the nucleons) coupled to 
fj and 7] in the presence of some external sources j (like e.g. photons or quark mass 
insertions) to the S-matrix (to be precise, we consider processes with one fermion line 
running through the pertinent Feynman graphs): 

1. Differentiation with respect to the external sources (here: j,r],f]). 

2. Multiplication with the physical propagator. 

3. Multiplication with the spinors u and u. 



4. Substitution of \fZu by Uphys and similarly for u. 


This leads to the form of eq.(3.10). Thus, Z]^ is dehned by the product of the physical 
propagator and the two-point Greens function. For our case, this leads to 


= (1 + Go-'Po)To-'n+ + (1 + 
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which means that the Z-factor appears as the sum of two terms, which are proportional 
to the velocity projection operators. The part of relevance for the present discussion 
is, of course, since all Greens functions corresponding to the positive velocity 
eigenstates include a projector P+. This is exactly what has been used in eqs.(3.9,3.25). 


B. Non—renormalization of the electric charge due to loops 


\ 

I 

/1 


a) 


b) 


c) 


\ 

jyv/vAz 


d) 


e) 


Figure 2: a) Self-energy graph. 
The solid and the dashed line 
denote the nucleon and the pion, 
in order. b)-e) One loop graphs 
contributing to the isovector 
charge form factor. The wiggly 
line denotes the photon. 


In this appendix, we explicitely show that the loops 
do indeed not renormalize the charge, which fol¬ 
lows from gauge invariance, but only lead to a mo¬ 
mentum dependence of the electric Sachs form fac¬ 
tor. Here, we entertain the possibility of a differ¬ 
ent prescription how to treat the spin vector in d 
dimensions. To be precise, as in ref. [4] we con¬ 
sider the Dirac form factor since at zero momen¬ 
tum transfer, P('(0) = ^^(0) because of the rela¬ 
tion GE{k‘^) = Fiik"^) — {k"^/4:m)F2{k‘^). However, 
the calculation shown it what follows differs from 
the published version [4] in that the spin matrices 
Sp are not extended to d 7 ^ 4 dimension, but only 
the loop momenta /. Such a procedure can also be 
extended to higher orders, if needed. For Sp, we 
consistently use the Pauli matrices, Sp = (d,d)/2. 
Therefore, a ■ I picks out three of the d — 1 space¬ 
like components of 1. Of course, the method used in 
refs. [4] [13] is correct, we only show this alternative 
way to demonstrate that while the Z-factor depends 


on such choices, physics does not (as long as one cal¬ 
culates correctly). Glearly, the loop contribution to the Z-factor in this approach is 
different to the one given in section 4. Gonsider hrst the self-energy graph in £g.2a. 
It gives (note that to this order we can identify all quantities in the chiral limit with 



their physical values, the error being of higher order) 

^ Ssi I dH .9g\ ,, , 

i“‘'( ) Fi i {2i,Y{Ml-r-){vl-Lj) 

using = —3/4, S' ■ u = 0, cu = u ■ /c, and the loop function J 2 {uj) is given in app. B 
of [13]. This leads to the well-known mass shift (which we already derived in section 3, 
but it is instructive to show that this somewhat unusual treatment of the spin matrices 
indeed leads to the correct result) 
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and the momentum-independent loop contribution to the Z-factor 


z'r = 1+su(o) = 


MMi Mi 

327r2F2 A 


modL(A) , 


(B.2) 


(B.3) 


where we are not concerned with the inhnite piece ~ L in what follows since it cancels 
as in described in [4] and the prime denotes differentiation with respect to cn. To arrive 
at this result, we have used ^ 2 ( 0 ) = ~^n = ln(M 7 r/A) (dropping again the 

term ~ L). Note the difference to in [4], where one has (31n(M7r/A) + 1) instead 
of 31n(M7r/A) (in a short-hand notation). In particular, in this formulation only the 
logarithmic terms survive. We do not absorb these in the LECs as done before. We 
now consider the charge form factor at = 0. The graphs lb,c contain for 

completeness we give the result for the isoscalar as well as the isovector form factor 


Ampdft + k) = - 1) = -(1 + r,)^^ 1.1 ^ 


The two one-loop diagrams Id and le give 
Amp(ld) = J 2 




Amp(le) = -Tj ^ ps(0) = n In — , 

using 78 ( 0 ) = (see app. B of [13]). Adding np pieces, we end np with 


(B.4) 


(B.5) 


(B.6) 


Amp(16 + Ic + 111 + le) ^ In ^ 


—3 — 3t3 -|- 3 — T 3 -|- 4t3 — 0 , (13.7) 


which is nothing bnt the anticipated resnlt. The terms ~ L, which we did not give, also 
cancel. But this is a different statement, since the cancellation of the inhnities only 
means that one has renormalized properly, and thus it does not constrain the hnite 
pieces we have discussed here. Finally, we remark that the considerations presented in 
this section can easily be extended to the next order. For the sake of brevity, we do 
not spell out the pertinent details here. 
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